Duality in the Azbel-Hofstadter problem and the two-dimensional d-wave superconductivity with a 

magnetic field 
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A single-parameter family of lattice-fermion model is con- 
structed. It is a deformation of the Azbel-Hofstadter problem 
by a parameter h — A/t (quantum parameter). A topological 
number is attached to each energy band. A duality between 
the classical limit {h — -\-0) and the quantum limit {h — 1) is 
revealed in the energy spectrum and the topological number. 

The model has a close relation to the two-dimensional d- 
wave superconductivity with a magnetic field. Making use of 
the duality and a topological argument, we shed light on how 
the quasiparticles with a magnetic field behave especially in 
the quantum limit. 



Two-dimensional Dirac fermions with a gauge field are 
of current interest, e.g., in the context of the vortex state 
in a two-dimensional d-wave superconductivity In 
our study, Dirac fermions with a gauge field are realized 
on a two-dimensional lattice. It is a single-parameter de- 
formation of the Azbel-Hofstadter problem [^|J|] . In this 
paper, the parameter h is called quantum parameter 
A topological number is assigned for each energy band 
10-^. As the quantum parameter is varied continuously 
from the classical limit (h = -1-0) to the quantum limit 
{h = 1), the energy spectrum is reconstructed through 
the change of each topological number (plateau transi- 
tion [^-11 1). Although the two limits are not connected 
adiabatically, we found that there is a duality between 
the classical and the quantum regime. 

The model has a close relation to the two-dimensional 
d-wave superconductivity with a magnetic field. Apply- 
ing the duality and a topological argument, we provide 
insights into the quasiparticle spectrum. In the quantum 
limit, interference effects become relevant especially at 
zero energy j|] . The existence of edge states is discussed 
as well. It reflects a non-trivial topology of each energy 
bandgll. 

Let us define a key Hamiltonian in our paper, which 
is a single-parameter family of lattice fermion model. It 
is a deformation of the Azbel-Hofstadter problem. The 
Hamiltonian H = J2i m c|W/mC„i with 
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Ai,n = -A„i/eR and 



Aim/^TT = (j) = pI q [p and q are coprime integers and 



the summation runs over four links around a plaquctte). 



Here LmeZ\ f' 



t. 



-AO 

^m±(0,l),m 



"m±(l,0),m 



''m±(l,0),m ~ ''m±(0,l),m 

A (<, AeR) and the other matrix ele- 
ments are zero, t is set to be a unit energy and a relevant 
parameter, quantum parameter, is defined by ft, = A/i. 
The relation of this model to the two-dimensional d-wave 
superconductivity with a magnetic field is discussed later. 
In the classical limit h = +0, this model decouples to 
two essentially equivalent Hamiltonians. It is the Azbel- 



Hofstadter Hamiltonian J2i m'^i^ '""^im^rn- The energy 
spectrum at = is given hy E = ±^y A{ky + \B{k)\'^ 
where A{k) — 2t{cos kr^ +cosky), B{k) = 2A(cosfca; — 
cosky) and tt, tt] x [— tt, tt]. The upper and lower 
bands touch at four points {±n/2, ±7r/2) in the Brillouin 
zone. The low-lying excitations around the gap-closing 
points are described by massless Dirac fermions. 

One of the basic observables is a topological number 
for the n-th band, C„ |@-§. It is 

Cn y"dk z • (V X A„), A„ = (w„(k)|V|w„(k)) 

where V = d/dk and |it„(k)) is a Bloch vector of the n-th 
band. The integration J dk runs over the Brillouin zone 
(torus). The non-zero topological number results in the 
existence of edge states. In order to see it, put the sys- 
tem on a cylinder and introduce a fictitious flux through 
the cylinder (it is equivalent to a twist in the boundary 
condition) p3[ . The edge states move from one boundary 
to the other as the fictitious flux quanta hc/e is added 
adiabatically. The number of carried edge states coin- 
cides with the summation of topological numbers below 
the Fermi energy |^,^ . Due to the topological stability, 
a singularity necessarily occurs with the change of the 
topological number (plateau transition |p|-|rT|]). The sin- 
gularity is identified with an energy-gap closing on some 
points in the Brillouin zone. 

In Figs. 1-3, the energy spectra are shown. As h is 
varied continuously from the classical limit (h — +0) 
to the quantum limit {h = 1), the energy spectrum is 
reconstructed through the plateau transitions. Although 
the two limits are not connected adiabatically, a main 
feature in the data is that there is a symmetry between 
the classical and the quantum regime. It leads to a claim 
that 
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' there is a faithful correspondence between cj) = (1/2 + 
p/q) in the classical limit h = +0 and (p = p/q in the 
quantum limit h = V. 

We call this phenomena duality. As discussed above, 
this model reduces to a doubled Azbel-Hofstader problem 
in the classical limit h = +0. The duality is an analogue 
of a statistical transmutaion (composite fermion picture) 
in the fractional quantum Hall effect In the com- 

posite fermion picture, a locally attached flux to each 
fermion is replaced with a global uniform flux. In our 
case, pairing effects (off-diagonal order) play the role of 
shifting a flux globally. It is reminiscent of a symme- 
try between the d-wave pairing and the 7r-flux, in other 
words, the SU{2) symmetry [|5|. Here some comments 
are in order. At </) = and 1/2, we proved analytically 
that the claim is exact. Moreover, as shown in Figs.l and 
3, topological numbers are consistent with the claim. 

As an application of the duality, let us study the two- 
dimensional d-wave pairing model with a magnetic field 
especially in the quantum limit {h = 1). The pairing 
model is H ^ 
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Under the unitary transformation c„|— >d„|, c„|^rf|^ 
(for Vn), it becomes H = Z]i,mMlT*'™'^™T + Ai'^im.dm.i + 
d\^Aimd^^^ -t- dmiA*i^di^]. It is equivalent to the 
Bogoliubov-de Gennes Hamiltonian for the singlet super- 
conductivity. Here t^^^ — tmi (hermiticity) and Aim — 
Ami (SU(2) symmetry) are imposed as well. It satis- 
fies a relation —{dyHimCy)* = Him, which results in a 
particle-hole symmetry in the energy spectrum. The two- 
dimensional d-wave pairing model with a magnetic field 
is defined by the pairing model (||) with 

tim = eyip{iAim)t°m^ 

Aim = ex.p{iipim)A°m 



In the following, we set a period for Aim and vim to 
be IxXly Here we emphasize that, in the context of 
superconductivity. Aim and vim are determined in a self- 
consistent way. Moreover, the spatial variation of |A;m| 
plays a crucial role especially near a vortex core. In our 
study, however. Aim and vim are treated as adjustable 
parameters and focus is put on the duality or topologi- 
cal arguments which do not depend on the detail of the 
potential. 

In Fig. 4, an example of the density of states (DoS) 
is shown for the d-wave pairing model with a magnetic 
field in the quantum limit (h — 1). Although the weak- 
field regime (0^0) is relevant to reality, we show the case 
0=1/5 for clarity. In the weak- field regime, the number 
of energy bands increases but the following arguments 
are robust. 

As the system at = p/q approaches the quantum 
regime (h^l), the quasiparticle spectrum becomes close 
to that of the Azbel-Hofstadter problem at (/) = (1/2 -|- 
p/q). The Landau levels are mixed due to lattice effects 
and Vim- It causes a singularity at zero energy |^ and 
the broadening of each level ('Landau bands'). It is due 
to quantum interference effects through spatially varying 
potentials. It is an analogue of the vanishing DoS at 
zero energy in random Dirac fermions [|l^-|l8|. In the 
case oi (j) — p/q (q=odd), it is a natural consequence 
of the duality. In other words, it can be interpreted by 
the fact that there exists a singularity at zero energy 
i.e. p{E = 0) = in the Azbel-Hofstadter problem at 
(f) — {\/2 + p/ q). It is also to be noted that, apart from 
the singularity and the broadening, the energy spectrum 
of the Azbel-Hofstadter problem at (j) = {1/2+p/q) {p = 
1, g^l) is '^±ujH\/n near the band center and ^zt{u!n + 
C) near the band edge (neN). 

Now put the system on a cylinder (periodic in the y 
direction and open in the x direction) and let us study 
the basic properties of edge states in the d-wave pairing 
model with a magnetic field. The Schrddinger equation 
for the Hamiltonian (0) is 



where Aim = -Ami, Vim = (</?; +¥'m)/2eR so that t,*^ 
tml (hermiticity) and Aim — Ami (SU(2) symmetry) are 
satisfied ||l|,|| . Performing a gauge transformation 



we obtain 
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where vim = (v/ — 'Vm) /2 — Aim- It is a lattice realization 
of the Hamiltonian discussed in ref. ||l| and vim corre- 
sponds to the superfluid velocity. In the case vim = 0, 
it reduces to the Hamiltonian dl^) and the quasiparticle 
spectrum consists of 'Landau levels' E = ±ujH\/n (nGN) 
in the continuum limit [III. 
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Decompose all the sites Af into two sublattices A and B 
{ALiB = Af and AnB = cj)) where tij and A^ connecting 
the same sublattice are zero. Define Wk by Wk = +Wk 
(fcsA) and Wk = —Wk (kGB). Then it follows from the 
SU(2) symmetry that 
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-f* 



= E 



It leads to the claim that, if an edge state exists, a paired 
(degenerate but linearly independent) edge state can be 
constructed which is localized in the same side |l9| . 

Next we shall discuss the existence of edge states. At 
first, set (j) = p/q. Vim = and focus on an energy gap 
with a non-zero topological number (see, for example. 
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Fig. 3. The duality implies that a topological number for 
a generic gap is non-zero even number in the quantum 
regime. It is consistent with the previous result j^]). As 
discussed above, a non-zero topological number leads to 
the existence of edge states in the gap. In other words, 
the existence of edge states is topologically stable. Next 
consider the case when vim is finite. As the vim is var- 
ied from zero to the finite value, the Landau bands are 
deformed and can be overlapped like a semimetal. We 
can, however, observe edge states due to the topologi- 
cal stability, when the plateau transition is absent in the 
deformation. In fact, see Fig. 5. The existence of edge 
states can be confirmed in each energy gap. It reflects a 
non-trivial topology of each energy band. In the above 
discussion, we have employed a cylinder. It is possible to 
consider the same problem on a geometry with 'defects' 
(e.g. annulus). In the case, the states analogous to the 
edge states on a cylinder may be bound to the defects. 

In summary, a single-parameter family of lattice- 
fermion model is constructed in two dimensions. The 
parameter h = A/t is called quantum parameter. A du- 
ality is revealed in the model between the classical limit 
(h = +0) and the quantum limit {h = 1). Employing the 
duality and a topological argument, we provide insights 
into how the quasiparticles with a magnetic field behave 
especially in the quantum limit. A more detailed study 
including a self-consistent potential and the fluctuation 
is left as a future problem. It is crucial for the analysis 
of dynamical properties. As discussed in the context of 
the integer quantum Hall effect, when the potential is 
sufficiently strong, it can bring the system to a totally 



different state |20 
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FIG. 1. (j) — p/127, h = A/t — 0. An integer in an energy 
gap is the summation of topological numbers below the gap. 



FIG. 2. (f) = p/127, h = A/t = 0.5. An integer in an 
energy gap is the summation of topological numbers below 
the gap. 



FIG. 5. Energy spectrum of the d-wave pairing model 
(h = A/t = 1) with a magnetic field {p = 1, q = 3, (p = 1/3) 
on a cylinder (the boundary is periodic in the y direction and 
open in the x direction). The spatial average and variance of 
Vim is and 0.1 respectively. The period IxXly is 30x3. The 
system is composed of 1x100 cells. 
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FIG. 3. (f) ~ p/127, h = A/t — 1. An integer in an energy 
gap is the summation of topological numbers below the gap. 



FIG. 4. An example of the DoS for the d-wave pairing 
model {h — A/t = 1) with a magnetic field {p = 1, q = 5, 
<j> = 1/5). The spatial average and variance of vim is and 
0.1 respectively. The period IxXly is 5x5. The DoS for the 
d-wave pairing model {h — 1) without a magnetic field is also 
shown (broken line). 
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